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Abstract 



In this work we develop an approach to obtain analytical expres- 
sions for potentials in an impenetrable box. It is illustrated through 
the particular cases of the harmonic oscillator and the Coulomb poten- 
^ ' tial. In this kind of system the energy expression respect the correct 

■ quantum limits, which is a very important quality. The similarity 

of this kind of problem with the quasi exactly solvable potentials is 
explored in order to accomplish our goals. 

Quantum systems under non trivial boundary conditions, corresponding 
to penetrable and nonpenetrable walls, do simulate the effect of atoms or 
molecules in the neighbor of a central particle. The dependence of the eigen- 
values on the box size, allows one to define the effect of the pressure over 
the system. Some examples would include that of the proton-deuteron trans- 
formation as a source of energy in dense stars jy], the determination of the 
escape rate from galactic and globular clusters 0, and the understanding 
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of the spectral line shift under pressure pi. One good list of phenomena 
associated to this kind of physical system is presented in [[|. Recently, the 
fabrication of the quantum dots in the semiconductor physics did brings a 
renewed interested on this kind of problem ||. Another example of applica- 
tion of this type of system, would be that of the Casimir effect and even that 
of Rydberg atoms between metallic plates . 

These important systems are treated by imposing the vanishing of the 
wave function in some finite point, once that the confined particle must 
be found in a finite region of the space. On the other hand, a numerical 
treatment of this problem has the natural drawback of becoming difficult 
any qualitative analysis, and even estimation, due to the need of calculating 
the energy for each new value of the potential parameter |7|], ||. Another 
possibility of studying the system is to apply to the Rayleigh-Schroedinger 
perturbation theory M, [TOI, variational methods |TT|, 



or some semi- 



classical approximation |]T3f. Recently it was applied the strong-coupling 



expansion approach to this kind of problem ||14|| , and considering the case of 



the pressure exerted by a stack of membranes upon enclosing walls [[HJ, [IR 



it has been treated particularly the case of a single membrane [17|, |18 



Here we show that the imposition of the vanishing of the wave function 
at the boundaries imply that the potential can be treated in a similar fashion 
to that of the so called quasi exactly solvable potentials |19[-||21[, which are 
quantum potentials which have part of the energy spectrum exactly solvable, 
provided that some relations between the potential parameters hold. In fact, 
we use this feature in order to deal with the problem in the same way that it 
was done in the case of an anharmonic oscillator , [2]|] . As an application 
of the method, we present an analytically approximated expression for the 
quantum problem of one particle under the action of a harmonic oscillator 
bounded to a finite region of the space. Furthermore, our analytical expres- 
sion does have all the correct physical limits as, for instance, that for large 
values of the principal quantum number. On the other hand, in contrast with 
other approaches, our result is valid for arbitrary values of the box size. 

First of all, we illustrate the method by treating the harmonic oscillator 
bounded in a "one-dimensional box" of length a. Some similar problems were 
treated previously in the literature 0, [12|], but in general only solutions 



valid for a specific value of the potential parameter are taken and, every 
time one needs the energy spectra for another value of that parameter, the 
numerical calculation must be repeated. This implies that one hardly can 
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visualize a qualitative behavior of the system, besides there is the obvious 
unpleasant need of repeated numerical calculations each time one needs to 
change the value of the potential parameter. By the other hand, methods like 
the perturbation theory are intrinsically limited due to the need of restricting 
the potential parameter or the length of the box. 



1 Harmonic Oscillator in a one dimensional 
box 

We start by treating the harmonic oscillator in a box, similarly to the treat- 
ment given to the quasi-exactly solvable potentials |19|- ||21|| , which did lead 



us to obtain analytically approximated solutions for the eigenvalues of an- 



harmonic oscillators under usual boundary conditions \22 \. 

Remembering that the general solution for the harmonic oscillator is given 
by something like 

(x) = P n (x) exp (-ax 2 ) , (1) 

where P„(x), in principle, is a polynomial of infinite degree (the parabolic 
cylindrical functions), which however is taken finite (Hermite polynomials) 
in the case of boundary conditions at the infinity, in order to get a conver- 
gent function. Here however, the exponential does not guarantees that the 
wavefunction vanishes at the walls of the interval, so that this condition de- 
termines the energy spectrum. Notwithstanding, when we do insist to using 
finite polynomials, which are nothing but the known Hermite ones, whose en- 



huj (n + |) , we get an equation for 
ved, give us those frequencies which 



ergy is well defined and given by E n 
the oscillator frequency which, once so 
generate polynomials whose zeroes are at the walls of the box. At this point 
it is important to stress the similarity of this problem with that of the so 
called quasi-exactly solvable potentials, in the sense that in this case only for 
some specific frequencies the eigenfunction will be elementary or, in other 
words, it is represented by a finite polynomial, the Hermite ones, times an 
usual exponentially decaying factor. Once one gets those frequencies, one 
must to substitute them on the above expression for the energy, so obtaining 
the corresponding exact energy for the harmonic oscillator in a box. 
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In order to obtain the frequency and the energy for an oscillator in a box 
whose walls are at -a/2 and a/2, by applying the idea which was outlined 
above, we use the condition 



H n (±0 = 0, (2) 

so that H (x) = 1, does not give us any solution, the polynomial H 2 (if) 

generates one solution, if 4 f±§) two solutions, so for and so on. To be more 

precise, let us exemplify with case of H 2 (±f). In this case the equation 
which implies that the wavefunction vanishes at the boundary is given by 

H 2 (±£) = M-2 + 40 = 0, (3) 
where £ = \J^jf-- Remembering that in this case E 2 = ^hcu, one gets finally: 

2H TP K ^ (AS 

ma z ma z 

As can be easily verified, the wavefunction of the harmonic oscillator for 
the above frequency and energy vanishes at the boundaries, not presenting 
any nodes between the walls. So it corresponds to the ground state of the 
oscillator in a box with impenetrable walls. Now, if one takes the following 
even case of (±§)j it is straightforward to show that there are now two 
possible solutions of ||, respectively corresponding to the ground state and the 
second excited one, so for and so on. Analogously a similar structure appears 
in the odd cases. The first solution of each polynomial always corresponds 
to the first excited state, the second to the third excited state, etc... 

However, before to do such a calculation, we take profit of a scaling 
symmetry of this kind of problem in order to simplify it. In doing so, we 
show that this problem has only one independent variable for the energy 



dependence. By choosing x = | y^j-y, and redefining the frequency and 
the energy conveniently, one gets 

where the new variables in terms of the old ones look like 

m a 2 ma 2 
u a = --u,E a = J¥ E, (6) 
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and now we need to find the zeros of ip such that ip (±1) = 0. The resulting 
solution will be valid for arbitrary values of the box length and its frequency. 
In this way, by solving a sufficiently great number of polynomials, we can 
plot the energy as a function of the frequency for each one of the energy 
levels, and then try to get a functional relation between these several levels. 
So obtaining an approximate analytical expression for the dependence on the 
frequency, the length of the box and the quantum number. 

At this point, in order to help our quest for a suitable analytical expression 
for the energy, we remember that in the limit of zero frequency one recalls the 
free particle in a box. On the other hand, when the length of the box goes 
to infinity, one should recover the well known harmonic oscillator spectra. 
Finally for great values of the principal quantum number, once again the 
energy becomes close to that of a free particle in a box. As a consequence of 
these physical constraints, we write the expression for the energy as 

8 = E = n2 + kuJa + gM ' (7) 



where 



1 

7T 2 V" 2 

and the function g (cu a ) should have the following appropriate limits 



k = ~2 ( n - «),"■= 1,2,3,... (8) 



g ( Wo = 0) = 0, g {u a -> oo) = 0. (9) 

Furthermore, one should also have this function going to zero when n in- 
creases or at most growing more slowly than n 2 . Using this hypothesis as a 
guide we try the following form for g n (u a ): 

g n {u a ) = c Q {n)u a e-^U^K (10) 

After performing the necessary calculations we verify that indeed, this 
function has the expected behavior supposed above. Besides, by obtaining 
the roots of Hermite polynomials up to two hundred degree, we did find 
an approximated expression for the coefficients Cj (n), taking J = 3, whose 
expressions are given below: 

Co(n) = (0.405231 + 0.810579 n)' (H) 
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, , 0.0104832 0.00588616 0.00187449 , , 

c i M = Z 3 3 > (!2) 



n n 2 n 3 



c 2 ( n ) = 10~ 6 [-1.24 + 1.35 cos h (0.1762 n - 0.12)]" 1 , n > 1, (13) 

c 3 ( n ) = 10~ 8 [-2.5 + 2.6 cos /i (0.086 n - 0.278)] -1 , n > 3, (14) 

Note that for the last two coefficients, the first elements were separated in 
order to getting better fittings. In these cases one have 

c 2 (l) = 3.70973 10~ 6 , 

c 3 (l) = -1.54146 10~ 6 , c 3 (2) = -8.78283 10~ 7 , c 3 (3) = 4.56951 KT 8 . 

(15) 

Note that, in fact, the coefficients really approaches to zero for larger 
values of n. These calculations were performed for n < 20, so that for higher 
values of n, one must extrapolate it, but it should be expected good results 
due to the behavior of g n (cu a ) when n becomes greater and greater. 

The comparison of the energy coming from the above analytical approx- 
imation with pure numerical values, in the case of the range of frequencies 
considered, shows that the error was always less than 0.07% along the range 
of the parameters verified which due to technique reasons, depend of the en- 
ergy level studied and grows for higher quantum numbers, below we present 
a Table with the corresponding ranges and their respective maximum per- 
centual errors for the first twenty levels. 



n 


W 1max 


/Omax 


n 




O /Omax 


1 


2.0 


9.0 X 10~ 5 


11 


60.5 


0.018 


2 


6.0 


0.072 


12 


67.3 


0.018 


3 


10.9 


9.1 x 10~ 3 


13 


74.1 


0.017 


4 


16.3 


0.02 


14 


81 


0.017 


5 


22.1 


0.016 


15 


87.9 


0.016 


6 


28.1 


9.3 x 10~ 5 


16 


94.9 


0.016 


7 


34.4 


0.012 


17 


101.9 


0.015 


8 


40.8 


0.017 


18 


109 


0.013 


9 


47.2 


0.019 


19 


116.1 


0.016 


10 


53.8 


0.019 


20 


123.2 


8.9 x 10~ 3 



It is interesting to stress that the type of energy expression obtained 
here, can be used to study the appearing of quantum chaos in the cases of 
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anharmonic systems [51], | 53| . This happens because one can use the energy 
expressions coming from this approach in order to get the number of levels 
for a given range of energy and the corresponding related parameters, used 
to discuss the appearing of the quantum chaos. 



2 Coulomb Potential in a sphere 



In this section we work with the more mathematically involved case of the 
Coulomb potential in a sphere of radius a. Notwithstanding, the basic pro- 
gram to be followed is the same of the previous section. Once again the 
general solution for the radial wavefunction of this problem is of the kind 

R(x) = K{x)e W {^-) (16) 

where V(r) = -a/r, and n is the principal quantum number. This polynomial 
will be of infinity degree, unless that we impose that the wavefunction van- 
ishes at the infinity, in which case the polynomial degree is finite, one gets 
the generalized Laguerre polynomials. However, as in the previous section, 
despite we are treating a case with finite boundary conditions for which it 
has in general an infinity number of terms, we will insist in using a finite 
polynomial in order to generate simple exact solutions, remembering that 
the price to be paid is that consequently, the energy is established and given 
by E n = —[lot 2 1 (2hn). As a consequence the constant a will be chosen in a 
manner to guarantee the vanishing of the wavefunction at the sphere surface. 
After that we substitute the so obtained constant in the above expression for 
the energy, getting the corresponding energy for the system. 

Starting from the radial Schroedinger equation, once that the angular 
part of it is unchanged due to the fact that the boundary condition preserves 



the spherical symmetry p9| : 

-h 2 1 d 2 



2/j, r dr 



[r(R(r))\ 



-\^ n +vir) 



R (r) = E R(r) (17) 



with V(r) = — — . Defining R(r) = we have 
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-h 2 d 2 l(l + i), 2 



2/1 dr 2 ^ 2pr 



k 2 + V (r) 



u (r) = E u (r) . 



(18) 



Redefining the radial variable as p = ^-r and then substituting it in the 
previous equation one gets 



C 



d 2 1(1 + 1) 



dp 



P A 



-E 



v(p)=0 



(19) 



4 

where C = By choosing the ansatz u(p) = v (p) exp (—ar) one obtains 



C 



d 2 



d 



, + 2a— - a z + 
dp z dp 



1(1 + 1) 2 



-E 



v (p) = 0. 



(20) 



P P\ 

By choosing v (p) = p l+1 J2o b n p n and performing the identification E = 
—Ca 2 , we get the following recurrence relation 

b - 2(l-q(n + / + l)) 

n+1 /(/ + l)-(n + / + 2)(n + / + l)° n - 

Taking into account the boundary conditions at the infinity, we can truncate 
the series at bjy (N, integer), getting b^+i = 0, so obtaining 

1 



N + l + V 



and as E = —Ca 2 , one gets the usual expression for the eigenenergies in such 
boundary conditions 



C 



pe 2 1 



En ~ ~(N+\+if - (22) 

We can still use the recurrence relation in order to obtain an expression 
for the generic coefficient b n in terms of bo. In doing so we get 



b ° 3=1 



' 2[a(j + l)-l] 
3(3 + 21 + 1) 



(23) 



As in the case of boundary condition at a finite radius, the series does not 
truncate and so we must to work with the infinite polynomial generated from 
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the above coefficients. This polynomial will be analogous to the Parabolic 
Cylinder function which appeared in the case of the harmonic oscillator, as 
treated in the last section. 

On the other hand the equation for R(r), can be written in terms of the 
Laguerre polynomial 



R(r) = r l exp (-At) L^{_ x (2kr) , (24) 



where k = \/^3r- Imposing that R(r = r ) = 0, we get 



-*&*®-0, (25) 

with g = 2/ ^ 2 r ° • In absolute analogy with the case of the harmonic oscillator, 
in this case we must to substitute the values of g coming from the above 
equation in that of the energy of the Coulomb potential, 

2ft V ~ Sn^rl [Zb) 



Now, by using the recurrence relation and imposing that the polynomial 
be finite, one obtains as expected the usual generalized Laguerre polynomials. 
In this case, the first solution coming from the equation |25] for each poly- 
nomial will give us the ground state energy, each second solution generates 
the first excited state, so for so on. Consequently we used the Mathematica 
software in order to perform this calculation in a systematic way. This was 
done for the first nine energy levels, using one hundred degree polynomials. 
For each energy level, we plotted the exact energies coming from the solu- 
tions described above as a function of the parameter g, and fitted it with a 
polynomial of the third degree given by 

E n ,i = E c^ l) g m (27) 

m=0 

where the coefficients C^' 1 ^ are presented in the Table 1. 
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n,l\m 



12 3 



1,0 
2,0 
2,1 
3,0 



-1.00421 0.502846 0.0152158 0.00487229 

-4.50591 0.918442 0.0402514 0.00162719 

-2.05892 0.384483 0.00167988 0.000448289 

-9.69984 0.938899 -0.0182804 0.000282908 



The energy expression coming from the use of the above parameters is 
in good accordance with the exact numerical data and, as happened in the 
harmonic oscillator case, does have a increasing range of validity with the 
increasing of the principal quantum number, as can be seen in the Table 
below. In these ranges, the error is quite small, about 10~ 2 %. 



By analyzing the data, one can verify that the range was limited by the 
maximum polynomial used (in this case we used polynomials of degree 100). 
So, in order to increase the precision and the range of validity, one should 
use higher degrees for the Laguerre polynomials. 

One interesting feature observed was that, in this case the energy is no 
more degenerate in the angular momentum quantum number, in contrast 
with the usual boundary condition at the infinite. So indicating the broke of 
some symmetry of the system. As the spherical symmetry was preserved by 
the boundary condition used, we must to look for another broken symmetry 
to be responsible for this behavior. In fact, the symmetry related to the 
conservation of the Runge-Lenz vector is that one. This can be argued by 
remembering that classically is associated to fact that the eccentricity of 
the elliptical orbits does not alter their energy. However in a finite sphere 
some of those will not be allowed. From the point of view of a geometric 
transformation, one can use the fact that the commutator of the Runge-Lenz 
vector with the radial coordinate, will generate a change in this coordinate 
such as 




gmin 0/0 max 



1.8 4.7xl0~ 5 

6.15 4.3xl0~ 3 

5.0 6.2xl0~ 4 

12.9 9.4xl0~ 3 



r' = r + {r, a.E}, {r, a.E} 



1 (rxL) 



a, 



(28) 



km r 
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where E and a are respectively the Runge-Lenz vector and the infinitesimal 
vectorial parameter of this symmetry transformation. From above one can 
conclude that when one is close to the frontier, the transformation could led 
it to beyond it, and this is not allowed by the boundary condition. Conse- 
quently this symmetry is now broken at the quantum level. It is interesting 
to observe that technically the Runge-Lenz operator and the Hamiltonian of 
the system is still commuting but, now due to the boundary conditions, the 
eigenfunctions of the Hamiltonian of this complete set of commuting opera- 
tors is not shared with those of the Runge-Lenz operator. It is interesting to 
note that this is a generic property of this kind of system. 

On the other hand it is not difficult to extend this approach in order to 
take care of other situations like, for instance, the case of finite walls, symmet- 
ric and antisymmetric boundary conditions, which can represent other kind 
of physical systems. Finally it should be interesting to study the applica- 
bility of this approach to systems with time-dependent boundary conditions 
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